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Abstract 

We review some recent techniques for dealing with non-hermitian random matrix models based on generalized 
Green's functions. We introduce the diagrammatic methods in the hermitian case and generalize them to the 
non-hermitian case. The results are illustrated in terms of the eigenvalue distribution, eigenvector statistics and 
addition laws. 

PACS: 02.10.Sp, 05.40.-a,24.60-k 
Keywords: matrix theory, random noise 



1. Introduction 

Non-hermitian random matrix models (NHRMM) 
permeate a number of interesting quantum prob- 
lems [Q,^, such as: open chaotic scattering, dissi- 
pative quantum maps, neural networks, quantum 
dots, non-hermitian localization, directed chaos, 
turbulence, QCD. 

In this short review, based mainly on our recent 
works, we go over some novel techniques for 
analyzing a large class of non-hermitian matrix 
models with unitary randomness (complex ran- 
dom numbers). 

As a main tool for investigating the properties 
of NHRMM we choose generalized, matrix-valued 
Green's function. We start by presenting stan- 
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dard arguments from Green's function theory for 
hermitian random matrix models. Then, using 
this construction as a reference frame, we gener- 
alize the idea to the case of NHRMM. Finally, we 
unravel several non-trivial properties of general- 
ized Green's functions, in particular the link with 
eigenvector statistics and non-hermitian addition 
laws. 

2. Hermitian RMM and diagrammatic ex- 
pansion 

The fundamental problem in random matrix 
theories is to find the distribution of eigenvalues 
Xi, in the large N (size of the matrix H) limit, i.e 



(1) 



where averaging (...) is done over the ensemble 
of N X N random hermitian matrices generated 
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with probability 



(2) 



where V is usually a polynomial in H. 

It is usually convenient to introduce the 
Green's function 



Giz) 



1 / 1 



(3) 



The eigenvalue distribution is easy to reconstruct 
from the discontinuities of the Green's function, 
using the relation 



^ = pi - inS{X) . 



X + ie X 
Indeed, 

- -Im lim G(A + ie) ^ ^ (Tr S{X - H)) 

N 
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1(^^(A-A,))=P(A). (5) 



This can be also recast in the form 
-d,G{z) = p(z) 

TT 
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where the derivatives are understood in the 'dis- 
tributional' sense. 

There are several ways of calculating Green's 
functions for HRMM. A number of methods are 
based on first diagonalizing the random matrices 
H UAW, then integrating over U using uni- 
tary invariance and finally analyzing the resulting 
integral over the eigenvalues. We will employ here 
the diagrammatic method, based on analogies of 
RMM to 0-fO dimensional field theory (see e.g. 
^). Its advantage is that it also works when the 
unitary symmetry is broken. This happens often 
in applications e.g. when the matrices have some 
block structure or are a sum of a random and a 
fixed matrix. We would like to emphasize that 
despite being apparently a perturbative method, 
the structure of the planar graphs that contribute 
allow for a resummation of the whole perturbative 
series and give the full exact result in the planar 
limit. 

A starting point of the analysis is the expres- 
sion allowing for the reconstruction of the Green's 



function from all the moments (TiH' 
~ ^ ^ N 
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We will use the diagrammatic method to eval- 
uate efficiently the sum of the moments on the 
right hand side. We note, however, that in gen- 
eral this series expansion is convergent only in a 
neighborhood of z = cxi. So the results of the di- 
agrammatic calculation apply directly only there. 
For hcrmitian random matrices this does not pose 
a problem. Since the eigenvalues lie only on some 
intervals on the real axis, by (||) the Green's func- 
tion is a holomorphic function of z on the com- 
plex plane except for the cuts on the real axis. 
Therefore we can reconstruct the Green's func- 
tion everywhere by analytical continuation. This 
last step is in fact trivial to perform for specific 
ensembles. 

For illustration, we consider now the well- 
known case of a random hcrmitian ensemble with 
Gaussian distribution. We introduce a generating 
fimction with a matrix source J: 



'rTrH^+TrHJ 



Z{J) = j dHe-- 
All the moments follow from Z{J), 



(Tri/") = ^— Tr ( ^ ) Z(J) 
^ ' Z(0) \dJ ' ^ ' 



,7=0 
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and are straightforward to calculate, since in 
our case the partition function reads Z{J) = 
exp2^Trj2. 

However, instead of calculating all the mo- 
ments individually alike by the Wick theorem, we 
could draw them using "Feynman rules" derived 
from our generating function, and perform a re- 
summation of all relevant graphs. The propaga- 
tor reads 



d^zjj) 

dJb„dJ'} 



1 dJ^ZjJ) 1 



The 1/z in (Q) is represented by a horizontal 
straight line. We depict the "Feynman" rules in 
Fig. I 
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The diagrammatic expansion of Green's func- 
tion is visualized in Fig. ||. Each "propagator" 
brings a factor of 1/iV, and each loop a factor of 
N. Therefore only planar graphs survive the large 
TV limit. Introducing the self-energy E compris- 
ing the sum of all one-particle irreducible graphs 
(rainbow-like), the Green's function reads 



(11) 



In the large N limit the equation for the self 
energy S, follows from resumming the rainbow- 
like diagrams of Fig. ||. The resulting equation 
( "Schwinger-Dyson" equation of Fig. ^ encodes 
pictorially the structure of these graphs and reads 



S = G. 



(12) 



Equations (11) and ( |12| ) give immediately G{z — 
G) = 1, so the normalizable solution for the 
Green's function reads 



G{z) = \{z-^/z^A) 



(13) 



which, via the discontinuity (cut) leads to 
Wigner's semicircle for the distribution of the 
eigenvalues for hermitian random matrices 



p(A) = ^V4^. 

3. Non-hermitian diagrammatics 



(14) 



Before we present a generalization of the pre- 
ceeding methods to non-hermitian random ma- 
trices it is instructive to understand what goes 
wrong when we apply the standard series ex- 
pansion (|^) to the simplest non-hermitian en- 
semble — the Ginibre-Girko one 0,^, with non- 
hermitian matrices A4, and measure 



P{M)=e 



(15) 



It is easy to see that all moments vanish 
(trA^") — 0, for n > so the expansion (|^) 
gives the Green's function to be G{z) = 1/z (di- 
agrammatically this follows from the fact that 
the propagator (A^gA4^) vanishes and hence the 
self-energy S = 0). The true answer is, how- 
ever, different. Only for \z\ > 1 one has indeed 
G{z) = 1/z. For \z\ < 1 the Green's function is 



nonholomorphic and equals G{z) = z. The reason 
for the failure is that in the expression con- 
figuration by configuration, the resolvent displays 
poles that are scattered in the complex z-plane. 
In the large N limit, the poles accumulate in gen- 
eral on finite surfaces, over which the resolvent is 
no longer holomorphic. Therefore we cannot an- 
alytically continue to the most interesting region 
where dG/dz ^ on the nonholomorphic surface, 
with a finite eigenvalue distribution. 

This points the way that a successful gener- 
alization should work directly for small z and 
that some regularization is necessary. Moreover 
we would not like to loose the calculational flex- 
ibility of diagrammatic calculations and express 
the quantities in terms of some (generalized) mo- 
ments. 

Exploiting the analogy to two-dimensional elec- 
trostatics, the following method |p|-pO[| has been 
proposed. Let us define the "electrostatic poten- 
tial" 

F = ^TY\n[{z- M){z- M'^) + e^]. (16) 
Then 



lim — „ ^„' — ^ = lim — ( Tr- 



:6 dzdz rA-a TV \ " (|z-7W|2 +£2)2 

Y,&^^\z-K))^7:p{x,y) 



(17) 



represents Gauss law, where z = x + iy. The last 
equality involves diagonalizing the non-hcrmitian 
matrix A^, by a linear {non-unitary) transforma- 
tion M = LAL^^. F can then be expressed 
in terms of the eigenvalues Ai, and the matrix 
V — {L''L)~^. a short calculation gives finally 
the required equality to the complex Dirac delta 
Siz-X,). 

In the spirit of the electrostatic analogy we 
could define the Green's function G{z,z), as an 
"electric field" 



~ dz 7Ve™\ \z-M^)iz-M) +e2) 



(18) 



However, instead of working ab initio with such 
quantity, and in view of applying diagrammatic 
and free-random variable methods it is much 
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more convenient to proceed differently as we now 
discuss. 

Following 1^,0 we define the matrix-valued re- 
solvent throughn 




z — M ie 
ie z - 

Qqq Qqq 
Qqq Qqq 

where we introduce the block trace defined as 



(19) 



TrB 



A B 
C D 



2Nx2N 



Tr A Tr B 
Tr C Tr D 



(20) 



Then, by definition, the upper-right component 
Qqq, is equal to the Green's function (p^). 

The block approach has several advantages. 
Let us define 



Z = 



z le 
ie z 



H 



M 








(21) 



Then the generalized Green's function is given 
formally by the same definition like the usual 
Green's function G, 



«4 



Tr, 



1 



z-n 



(22) 



More importantly, also in this case the Green's 
function is completely determined by the knowl- 
edge of all matrix-valued moments 



Tn 



z-^Hz-^n. 



(23) 



This last observation allows for a diagrammati- 
cal interpretation. The Feynman rules are anal- 
ogous to the hermitian ones, only now one has 
to keep track of the block structure of the matri- 
ces. To demonstrate this, we consider again the 
case of complex Gaussian randomness (the Girko- 
Ginibre ensemble) 

In this case the double line propagators are 



N 



d c ' 



(25) 



s ligh tly different realization was later proposed in 
and pa. 



As previously, we introduce the self-energy E, 
(here of course matrix- valued) , in terms of which 
we get 



Q = {Z- S)- 



(26) 



The resummation of the rainbow diagrams for S 
is more subtle. Instead of the hermitian equation 
S = G • = (due to the vanishing of the first 
propagator in (psj)), the analogue of ( p^ reads 
now: 



E = 



^qq ^qq 
Qq 



Qqq ■ Qqq 
Qqq Qqq ' 



fqq 







(27) 



The diagonal entries are zero, due to the struc- 
ture of the propagators ( p5| ) . The two by two ma- 
trix equations (26-2^ completely determine the 



problem of finding the eigenvalue distribution for 
complex randomness. Inserting ( |2^ ) into (|2^) we 
get: 



Qqq 

Q 



Qqq 
Qqq 



99 ='99 



1 



|zp - QqqQq 



='99 



^99 

Z 



.(28) 



Looking at the off-diagonal equation we see that 
there are two solutions: first one with 5^ = 
QqqQqq — 0, sccoud with b ^ 0. Tlic first one 
is holomorphic, and a straightforward calculation 
gives 



G(z) = -. 
z 



(29) 



The second one is nonholomorhic, and leads, via 
Gauss law, to 



G{z,z) 



p{x,y) 



1 d 



G{z,z) 



1 



(30) 



Both solutions match at the boundary b — , 
which in this case reads zz = 1. In such a sim- 
ple way we recovered the results of Ginibre and 
Girko for the complex non-hermitian ensemble. 
The eigenvalues are uniformly distributed on the 
unit disk |zp < 1. 

This simple example illustrates more general 
properties of matrix valued function. Each com- 
ponent of the matrix carries important informa- 
tion about the stochastic properties of the sys- 
tem. There are always two solutions for Qqq, one 
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holomorphic, another non-holomorphic. The sec- 
ond one leads, via Gauss law, to the eigenvalue 
distribution. The shape of the "coastline" border- 
ing the complex eigenvalues is determined by the 
matching conditions for two solutions, i.e. is de- 
termined by the equation QqqQqq = 0. 1 The prod- 
uct of the off-diagonal Green's functions hides an- 
other important information about the random 
ensemble. Recently, we have proven that 
the correlator between left and right eigenvectors 
(introduced in Q) is given by 



0{z) 



{La\La){Ra\Ra)5{z - Xa) 



-^G G 

yqqyqq ■ 



(31) 



For the case of the Girko-Ginibre ensemble the 
product QqqQqq foUows immediately from the ma- 
trix equations (|2^p7|) and reproduces the result 
111 



0(z) = -(l-|z|^ 



4. Free Random Variables 



(32) 



One of the powerful calculational features of 
generalized Green's functions stems from direct 
applicability of Free Random Variable (FRV) cal- 
culus to this formalism. Before we demonstrate 
this, we recall briefly that the concept of addi- 
tion law for hermitian ensembles was introduced 
in the important work by Voiculescu [ p^ . In brief, 
Voiculescu proposed the additive transformation 
(R transformation) , which linearizes the convolu- 
tion of non-commutative matrices, alike the log- 
arithm of the Fourier transformation for the con- 
volution of arbitrary functions. This method is 
an important shortcut to obtain the equations for 
the Green's functions for a sum of matrices, start- 
ing from the knowledge of the Green's functions of 
individual ensembles of matrices. This formalism 
was reinterpreted diagrammatically by Zee [p^ , 
who introduced the concept of Blue's function. 

Let us consider the problem of finding the 
Green's function of a sum of two independent 



^The shape can be also inferred from associated hermitian 
models using conformal mappings, for details see B. 



(free ||l^) random matrices Hi and H2, provided 
we know the Green's functions of each of them 
separately. First, we note that the IPI self-energy 
E can be always expressed as a function of G itself 
and not of z as usually done in textbooks. For the 
Gaussian randomness, 2(6*) = G (see (|2|)). Sec- 
ond, we note that the graphs contributing to the 
self-energy T,i^2{G) split into two classes, belong- 
ing to Si(G) and £2(6'), due to the independence 
of probabilities P{Hi) and P{H2) in the large A'' 
(planar) limit. Therefore 



I]i+2(G) = Si(G)+S2(G). 



(33) 



Note that such a formula is not true if the ener- 
gies are expressed as functions of z. Voiculescu 
R transformation is nothing but R{G) = E[G]. 
The addition ( |33| ) reads, for an arbitrary com- 
plex u, Ri^2{u) = Ri{u) -\- i?2(u)- The R opera- 
tion forms an abeli an g roup. The Blue's function, 
introduced by Zee |16|| , is simply 



B{G) = S(G) + G" 



(34) 



Therefore, using the identity G{z) = l/{z — S), 
we see that the Blue's function is the functional 
inverse of the Green's function 



B[G{z)] ^z. 



(35) 



In the case of hermitian Gaussian random matri- 
ces, R{z) — z, and B{z) = z + l/z. 

The addition law for Blue's functions reads 



Bi+2{z) = Bi{z) + B2{z) - 



1 



(36) 



The algorithm of addition is now surprisingly 
simple: Knowing Gi and G2, we find (|3^) Bi 
and i?2- Then we find the sum B1+2 using (p6[), 
and finally, we get the answer G1+2, by reapply- 
ing (^sl). Note that the method treats on equal 
footing the Gaussian and non-Gaussian ensem- 
bles, provided that the measures Pi and P2 are 
independent (free). 

The freeness of the ensembles, when rephrased 
in the diagrammatical language, means that: 
(i) Propagators and vertices corresponding to 
"Feynman rules" for Hi and H2 ensembles are 
disjoint, i.e. there are no vertices which link si- 
multaneously fields of type Hi and i?2- 
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(ii) The propagators cannot cross each other. 
Large N hmit ("asymptotic freeness") ensures 
the second condition by banning non-planar di- 
agrams. 

The Free Random Variable calculus in the guise 
of Blue's functions works extremely well for her- 
mitian random matrix models (see ) . However 
if we were to apply it to the Ginibre-Girko ensem- 
ble obtained as a sum of a gaussian hermitian Tii 
and antihermitian i7i2 ensembles, we would run 
into the same difficulties which plagued the stan- 
dard diagrammatic method. Explicitly the Blue's 
functions would read|^ 



Bh{z) 



1 1 



B^n{z) 



1 1 



(37) 



Addition law gives then BQinibre{z) = 1/z, lead- 
ing to the result G{z) = 1/z the same as that 
obtained by blindly applying 'hermitian diagram- 
matics'. This shows that one has to generalize the 
addition law in order to apply it to NHRMM. 

We demonstrated in the first part of these 
notes, that by introducing matrix-valued Green's 
functions we managed to extend the diagram- 
matical analysis to the case of NHRMM. There- 
fore, by imposing the rules (i) and (ii) on non- 
hermitian diagrammatics we can define the free- 
ness condition for non-hermitian ensembles. This 
analogous structure allows us to extend now with- 
out much effort the addition formalism to the 
non-hermitian caseQ. 

The generalized Blue's function is now a 
matrix valued function of a 2 x 2 matrix variable 
defined by 



B{g) = Z 



(38) 



where Z and Q were defined in ( pl| ) (the e in the 
definition of Z can be here safely set to zero). 
This is equivalent to the definition in terms of 
the self-energy matrix E, defined by the general 
relation (^, 



BiG) 



(39) 



^The factor 1/2 reflects just the necessary rescaling of the 
width of the gaussian. 

^Alternatively we could start from the generalized mo- 
ments and apply directly the generalized addition 
theoremsfor "symmetric moments" of Voiculescu. 



The diagrammatic reasoning as before leads to 
the addition formula for the self-energies and con- 
sequently for the addition law for generalized 
Blue's functions [|,| 



z = BiiG) + 62(5) - g- 



(40) 



For illustration we come again to the case of 
Ginibre ensemble. In this case the generalized 
matrix Blue's functions are 



6- 



H 



B. 



a b 
b c 

a b 
b c 



a b 
b c 



a b 
b c 



If -a 
2 I b 



a b 
b c 



(41) 



.(42) 



The generalized Blue's function of the Ginibre en- 
semble is thus 



B, 



Ginibre 



a b 
b c 



b 
b 



a b 
b c 



(43) 



Solving the defining equation ( pSD is equivalent to 
( p^ ) and, as we saw, gives the full answer. 

The power of the addition law for NHRMM 
stems from the fact that it treats Gaussian 
and non- Gaussian randomness on the same foot- 
ing [Q. Also it is a versatile way of obtaining 
results for non-hermitian random matrix ensem- 
bles starting from some simple building blocks, 
without having to go through explicit construc- 
tions of (supersymmetric) sigma-models for the 
relevant ensembles. In conjunction with ( ^l|) this 
method provides, in addition, information on the 
eigenvector statistics 0{z). 

5. Conclusions 

In this short review we tried to emphasize the 
versatility of the generalized Green's function ap- 
proach and underline the connections between the 
hermitian and non-hermitian treatment. Several 
details and practical applications of the results 
presented here are included in already published 
papers 
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Figures 



z "a ' 



(Him) 



AT "b"d 



Figure 1. Large N "Feynman rules" for Gaussian HRMM 
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Figure 2. Diagrammatic expansion of Green's function (^) for Gaussian ensemble. The 2"*^ and 4"* 
graphs are "rainbow" graphs contributing to the self-energy S. 




Figure 3. Schwinger-Dyson equation. 



